Effects of boundaries in mesoscopic superconductors 
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Abstract 

A thin superconducting disk, with radius i? = 4^ and height H = S,, is studied in the presence 
of an apphed magnetic field parallel to its major axis. We study how the boundaries influence the 
decay of the order parameter near the edges for three-dimensional vortex states. 
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A few years ago the response of a mesoscopic superconducting disk to a magnetic field 
parallel to its major axis was measured [l| and more recently the detection of giant vortices 
was made thanks to new advances in small-tunnel-junction technology \'^. The small volume 
to area ratio of mesoscopic systems brings new and interesting physical properties such as 
the onset of spontaneous and persistent currents ^] and for this reason the treatment of 
boundaries must be carefully considered. In this work the three-dimensional Ginzburg- 
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FIG. 1: l^'l^ vs. the distance from the center of the disk is shown here for the case of zero and one 
vorticity. The dots correspond to j^'p values at mesh grid points. The three-dimensional figures 
correspond to isocontours taken at 20% of its maximum value. 
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FIG. 2: Magnetization (arbitrary units) vs. the applied field. Inset shows the Hc2 to i?c3 region. 
The magnetization is normalized by a negative arbitrary constant. 

Landau (GL) theory is solved in a parallelepiped cell that contains a smaller mesoscopic 
superconducting disk inside. Space is discretized and gauge invariance kept on a 61 by 61 
by 16 grid. The distance between two consecutive mesh points along any of the major axes 
is ^/5, ^ being the coherence length at some fixed temperature. Thus the disk occupies the 
center of a 12.0^ by 12.0.^ by 3.0.^ cell. The vortex state solutions are obtained by numerical 
minimization of the GL free energy in the cell through the method of simulated annealing. 
Coupling of the disk to the outside non-superconducting space is included and the order 
parameter converges to zero outside the disk at the end of the minimization procedure. 
This makes the present method somewhat different to those that just seek the solution 
inside, but not outside, the disk [sl. The present description is restricted to a hard type 
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FIG. 3: Free energy vs. the applied field is shown here, the inset shows the Hc2 to region. The 
free energy is normalized to -ff^/Svr. 

II superconductor whose free energy is simply given by, / = / ^ (t^^|(V — 27ri/0oA)'?/'p — 
t\iI)\'^ + llV'l'^)! expressed in reduced units such that the density is normahzed to one. The 
magnetization is determined from M = const j dv :x. x J, where J is the supercurrent. The 
shape of the mesoscopic superconductor enters directly into the free energy through the step- 
like function t{x), equal to one inside the sample, and zero outside. As shown in Ref. |J] this 
approach yields the de Gennes boundary condition. Here we make the r(x) function smooth 
with an exponential decay linking the two sides. The condition of a vanishing supercurrent 
component pointing outside the superconductor is only approximately enforced between the 
disk and the outside external space: r = exp(— |(x — R)/?7|^). We take here for the two 
adjustable parameters, rj = 0.8^ and = 8. Numerically, this smoothness is advantageous 
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since a steep r(x) function directly couples a mesh point in the disk to an outside point 
rendering a severe depletion of the order parameter inside. In the present approach the order 
parameter never reaches its maximum value of one, predicted for a bulk superconductor, 
regardless of the 1] and values because the coupling to the outside conspires to lower the 
superconducting density near the edge of the disk, as exemplified in figure [T] for the cases 
of none and one vortex inside the disk. The magnetization and free energy curves obtained 
with the present approach display a total of 13 lines as shown in figures [2] and [3l respectively. 
Each line corresponds to a distinct vorticity state. The crossing of these lines in case of the 
free energy defines the so-called matching fields Gj. The present results are in fair agreement 
to those found by Baelus [t] for an extremely thin disk using a two-dimensional approach 
for the superconducting density. 

We have theoretically studied the vortex states of a mesoscopic disk with finite thickness. 
Previous studies have also considered a mesoscopic disk but with a vanishing thickness^]. 
We conclude that the present minimization procedure of the Ginzburg-Landau free energy 
is able to obtain the vortex patterns of truly three-dimensional mesoscopic superconductors. 
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